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I.  INTRODUCTION 


The  problem  of  realizability  of  biquadratic  functions  as  driving 
point  immitance  has  been  studied  by  a  number  of  investigators  [1]  -  [11], 
with  realizability  conditions  given  for  various  special  cases  as  well  as 
for  the  general  case  of  biquadratic  functions.  It  has  been  proved  [10] 
that,  given  a  pair  of  poles,  the  location  of  the  zeros  are  limited  by 
two  curves,  one  of  which  is  the  locus  of  the  zeros  such  that  the  biquad¬ 
ratic  function  is  a  minimum  positive  real  function. 

In  Chapter  II  of  this  paper,  various  realizability  conditions  are 
reviewed  and  illustrated.  FORTRAN  programs  are  written  in  Chapter  III 
for  illustration  as  well  as  comparison  of  different  cases. 

Many  forms  of  realizations  of  positive-real  biquadratic  functions 
are  well  known.  Some  of  these  forms  are  utilized,  in  Chapter  IV,  to 
study  the  sensitivity  problem  related  to  the  realizability  of  the 
function.  The  method  of  graphical  representation  by  Chan  [13]  is  applied 
in  this  study.  A  computer  program  is  written  for  the  illustration  of 
the  variation  of  sensitivity  phasors  [13]  as  functions  of  element  values. 

Finally,  in  Chapter  V,  a  discussion  on  the  results  is  made  together 
with  some  suggestions  for  further  studies. 
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II.  REALIZABILITY  CONDITIONS 


A.  INTRODUCTION 

Consider  a  general  biquadratic  function 


or 


Z(s) 


2  2 
s  +  2a  s  +  CO 
z  z 


2  2 
s  +  2a  s  +  CO 
P  P 


(s  +  z*)(s  +  z") 

(s  +  p’) (s  +  p") 


z(jco) 


2 

-CO 

2 

-CO 


+  2a  ICO  +  CO 
-zr 

+  2a  ICO  +  CO 
P 


2 

2 

P 


(co^  -  co^)  + 

2  2 

(co  -  CO  )  +  ico2a 
P  P 


(1) 


(2) 


It  has  been  investigated  by  Chan  and  Phung  [10],  that  a  biquadratic 

function  Z(s)  as  given  in  equation  (1)  with  conjugate- comp lex  zeros  in 

the  open  half  plane  is  realizable  as  a  driving  point  immitance  with 

passive  one-port  if  and  only  if  the  zeros  are  located  either  on  the 

curve  C  or  C.  or  inside  the  area  delimited  by  curves  C  and  C.,  whose 
0  1  0  1 

polar  equations  are  given  respectively  by 


max  CO  =  CO  +  2a  cos  <f)  +2  cos  (J)  (a  cos 


z  p  p 


and 


Min  CO  =  CO  +  2a  cos  , 
z  p  p  z 


z  p 


a  cos  (J)  (a  cos  ^ 
p  ^z  p  z 


*  +  CO  )"] 

z  p'J 

*  +  0)  )"] 

z  p  J 


(3-a) 


(3-b) 


where  (p  denotes  the  compliment  of  the  argument  of  the  zero  z  ,  such 
z 

that  , 

a  =  CO  cos  (p 
z  z  z 
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B.  NECESSARY  AND  SUFFICIENT  CONDITIONS 

To  obtain  the  real  part  of  Z(jcD)  of  equation  (2),  multiply  and 
divide  this  equation  with  the  conjugate  of  the  denominator,  so  that 


Re 


2  2  2  2  2 
(m  -  a))(a)  -  oi))+4aa(jj 

je z 


^  ,  2 
-0)  )  +  4a  (1) 

P 


(5) 


Z(s)  is  positive  real  if  the  numerator  of  expression  (5)  is  nonnegative 
for  all  values  of  m,  that  is 


.T/  2^  4  2,  2  ^  2  ,  .^22 

N(a))=(i)  -  (d((ji)  +0)  -4aa)+a)a) 

zp  zp  zp 


(6) 


Now  we  have  the  following  two  conditions  [10]: 

Condition-1 

2 

N  (o)  )  has  complex  or  equal  positive  real  roots. 

Condition- 2 

2 

N  (o)  )  has  nonpositive  real  roots. 

2 

Condition-1  implies  the  discriminant  of  NCm  )  is  negative  or  zero 

i.  e. 

D  =  (o)^  +  (1)^  -  4a  a  “  4(i)^(u^  ^  0  (7) 

z  p  z  p  z  p 
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or  after  factoring  we  get 


2  ir  2 

(o)^ 


Z  p  Z  pi  I  z  p  z  p 


,]<«  > 


thus  condition-1  is  satisfied  if 


2  2 
(oj  -  0)  )  4:  4a  a  <  (o)  -hu  ) 

z  p  z  p  z  p 


(7-a) 


(7-b) 


The  right  hand  inequality  is  shown  as  strict  since  inequality  corresponds 
to  real  roots  which  would  in  this  case  be  nonpositive,  which  is  a 
violation  of  condition-1. 

Next,  condition-2  implies  the  following  two  conditions  simultaneously: 
Condition-2a 

D  >  0  . 

Condition- 2b 

2  2 

aj+(jj-4aa^  0 
z  p  z  p 


It  can  be  seen  that  both  conditions  2a  and  2b  are  satisfied  if 


4a  a  X  ((1)  +  (1)  ) 
z  p  z  p 


(8) 


With  equations  (7-b)  and  (8)  we  obtain  the  necessary  condition  for  positive 
realness  of  Z(s)  to  be 


(o)  -  (jj  )  <  4a  a 

^  z  p  z  p 


(9) 


Now  suppose  that  equation  (9)  is  satisfied  with  the  equality  sign. 

2 

Z(s)  is  then  a  minimum  positive- real  function  since  N((jj  )  has  a  double 
positive  root  oj^,  where 


(Jj.  =  (1)  0) 

1  z  p 


(10) 
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Next  suppose  equation  (9)  is  satisfied  but  not  necessary  as  an  equality 


then 


a 

z 


4a 

P 


(11) 


The  necessary  and  sufficient  condition  obtained  in  the  above  discuss¬ 
ion  can  be  conveniently  and  geometrically  interpreted. 

Suppose  that  the  zeros  z*  and  z”  are  conjungate  complex  with  negative 
real  parts  and  the  poles  p*  and  p"  are  similarily  located  in  the  open 
left  half  plane. 


Let 


a  =  oj  cos  6 
z  z 


(12-a) 


where  4)  denotes  the  compliment  of  the  arg\nnent  of  the  zero  z*. 


If  0)  ^  0)  ,  then  we  let  m  m  +  u 

z  p*  z  p  o 


(12-b) 


If  0)  ^  oj  ,  then  we  let  m  m  -  u.  (12-c) 

z  p  z  p  1 

substituting  equations  (12-a)  and  (12-b)  into  equation  (9),  and  after 
rearranging  we  get 

2  '  * 

u  -  4a  (cos  4>  )  u  -  4a  a  cos  4>  <0  (13-a) 

o  p  O  P  P  2 


and  substituting  equations  (12-a)  and  (12-c)  into  equation  (9)  we  get 


u.  +  4a  (cos  *  )  u.  -  4a  a  cos  ({>  <0 

1  p  '  ^z'  1  p  p  z 


(13-b) 


Both  polynomials  in  the  left  hand  sides  of  the  inequalities  (13-a)  and 
(13-b)  have  a  positive  real  root  and  a  negative  real  root,  so  that  these 


11 


two  conditions  will  be  satisfied  if  and  are  nonnegative  and  smaller 
than  or  equal  to  the  positive-real  root  of  the  left  hand  side  polynomials 
in  (13-a)  and  (13-b)  respectively. 

Thus  solving  for  the  positive  real  root  and  then  substituting  the 

I 

result  into  equations  (12-b)  and  (12-c)  gives,  for  each  value  of  (()  >^  90°, 

z 

the  maximum  and  minimum  values  of  for  which  the  necessary  and  suffi¬ 
cient  condition  for  the  positive  realness  of  Z(s)  is  satisfied. 
Alternatively,  we  substitute  equations  (12-b)  and  (12-c)  into  equation 
(9)  and  get 

2 

u  -  4a  a  0 

z  p 

and 

2 

u.  -  4a  a  0  . 

1  2  p 

Solving  these  two  equations  for  the  positive  real  root  and  substituting 

the  results  into  the  expression  for  03^,  gives  the  maximum  and  minimum 

values  of  03  for  which  the  positive  realness  of  z(s)  is  satisfied, 
z 

Thus  we  have  in  this  case  the  two  equations 

max  03  =  03  +  2  /a  a  (14-a) 

z  p  z  p 

and 

min  03  =  03  -  2  /a  a  .  (14-b) 

z  p  z  p 

Note  that  these  equations  will  give  the  same  results  as  equations  (13-a) 
and  (13-b)  when  poles  of  Z(s)  are  complex,  and  will  also  give  the 
realizable  area  of  zeros  of  Z(s)  when  the  poles  are  real. 

Thus  given  Z(s)  with  complex  conjugate  poles  in  the  left  half  plane, 
or  with  negative  real  poles,  Z(s)  will  be  realizable  if  the  zeros  of  Z(s) 
are  complex  and  lie  on  or  within  the  curves  of  equations  (14-a)  and  (14-b) 
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In  the  special  case  where  Z(s)  has  real  poles  and  complex  zeros, 
it  has  been  shown  that  z(s)  can  be  realized  by  a  five  element  series 
parallel  circuit  [11]. 

The  key  restriction  in  this  case  is  given  by 

p”  (P"  -  p’)  -  (z’  -  p’')(z"  -  p")  >  0  (15) 

where  p',  p",  z'  and  z”  are  poles  and  zeros  of 


Z(s)  =  (s  +  z')(s  +  z") 
'  ^  (s  +  p’)  (s  +  p") 


and  p"  >  p' . 

If  we  substitute  a's  and  m's  for  the  p's  and  z's  of  (15),  then  we 

get 

I  P-  ^ 

p  =  -a  +  /a  -  to 

P  P  P 


fl  T 

-a  -  /a  -  to 
P  P  P 


^  p.  r 

-a  +  /a  -  to 


/2  T 

-a  -  /a  =  to 
z  z  z 


we  also  obtain 


+  2o  -  0)^ 


2a  a 

z  p  z  p 


2  2  ^  . 

u)  -  to  >0 
z  p 


rearranging  terms,  we  have 


2  2 

to  +0)  <  2a  a 

z  p  z  p 


+  2 

z  p  p 


Now  since 


2  2  2  2  2 
to  +(o  =(to  —  to)  +  2to  to  , 
z  p  z  p  z  p 
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then  we  get 


2  ^2  2 
(o)  -0)  )  ^  2a  a  +  2a  /a  -  0)  -  2a)  w 

zp  zp  zp  p  zp 


Comparing  this  condition  to  that  for  the  general  case  in  equation  (9) , 
we  obtain 


2a  a  +  2a  /a  -  0)  <  4a  a 

z  p  z  p  p  z  p 


2  fl  2 

(o)  -  0)  )  <  2a  a  +  2a  /a  -  o) 

z  p  zp  zp  p 


>  0)  <  4a  a 

z  p  z  p 


so  that  this  is  indeed  more  restrictive  than  the  general  case  and  we 
should  not  be  surprised  to  find  the  area  of  realizability  for  this  case 
enclosed  by  that  of  the  general  case. 

In  fact  this  argument  could  be  used  as  a  proof  of  the  realizability 
of  Z(s)  given  in  (15) ,  since  the  area  of  this  case  is  always  enclosed 
by  that  of  the  general  case  as  shown  above. 

The  construction  of  the  realizability  region  in  this  case  will  be  in 
accordance  with  reference  [11] ,  that  is 


i  i 

pM-r-p  (p-p)  , 


where  p"  >  p’. 


C.  SAMPLE  ILLUSTRATIONS 


(a)  Complex  poles  >  a^) 


Let  a  =  3 
P 


and  0)  =  5  . 

P 


14 


Then  from  equations  (14-a)  and  (14-b)  we  get 


and 


(jj  max  =  m  +  2/a  a 
z  p  z  p 


=  5  +  2/3^"  =  5  +  /12a  , 

z  z 


(D  min  -  0)  -  2/a  a 

z  p  z  p 


5  -  /12a  . 

z 


We  obtain  the  limitation  of  a  as  follows: 

z 


a^  <  for  complex  zeros. 


If  we  let  0)  max  =  a  ,  then  we  get 
z  z 


u  =  5  +  /12a 
z  z 


=  20.8  =  a  ,  maximum  value  for  a  ,  for  t 
z’  z 


Similarly  if  min  =  a^,  then 


0)  =  1.2  =  a  ,  maximum  value  for  a  ,  for  m 

z  z  z  : 


(b)  Real  Poles 


Let 


a  =2.5  and 
P 


0)^  =  2  with  poles  at  (-1,  -4). 


Then  we  obtain 


(1)  max  =  0)  +  2/  a  a 

z  p  z  p 


I  max. 
z 


min. 


=  2  +  /10a 


and 


iii  min  *  0)  -  vlOa 

z  p  z 

If  0)  max  =  a  ,  then 
z  z* 

(1)  =  2  +  /lOa 

z  z 

=  13.7  =  a^,  max  for  o)^  max  and  complex  zeros. 

From  the  examples  given  in  this  section  it  can  be  seen  that  the 
realizability  region  of  a  biquadratic  function  is  dependent  upon  the 
restriction  we  place  on  the  function.  That  is,  as  we  place  more 
restrictions  on  the  function  and  proceed  from  a  general  function  with 
complex  poles  to  one  with  real  poles  which  is  realizable  by  certain 
technique,  we  see  that  at  each  step  the  realizability  region  decreases. 
This  decrease  of  freedom  from  the  additional  restriction  is  to  be 
expected. 
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III.  COMPUTER  STUDY  ON  REALIZABILITY 


For  a  given  location  of  the  complex  pole  pair  (a  and  w  specified) 

f 

in  the  open  left  half  plane,  and  are  simply  plots  of  versus 

Thus  the  computer  program  was  written  to  accept  and  as  inputs  and 

to  generate  and  on  either  the  line  printer  or  the  Calcomp  plotter. 

Some  sample  outputs  are  shown  in  Figure  A-l.A  and  A-l.B  in  Appendix  A. 

The  curves  are  for  a  normilized  set  of  pole  locations  with 

0)  =  1 
P 

and 

0°  <  ((.'  <  90°, 

P 

?  T 

where  <b  is  the  compliment  of  the  argument  of  the  pole  p  . 

P 

Varying  (f)^  from  0^  to  90°  corresponds  to  shifting  the  poles  along 
a  circular  path  from  the  negative  real  axis  to  the  jm  axis. 

The  program  was  then  modified  to  produce  a  set  of  curves  and 

for  various  pole  locations.  Given  a  specific  pair  of  poles,  two  types 
of  plots  may  be  obtained,  depending  upon  the  requirements  of  the  problem. 

If  an  accurate  plot  of  the  realizability  region  is  required,  the 
program  can  be  used  so  that  a  plot  of  the  realizability  region  is 
required,  the  program  can  be  used  so  that  a  plot  similar  to  that  shown 
in  Figure  A-l.A  is  obtained.  Otherwise  if  a  quick  estimate  is  all 
that  is  desired,  a  plot  similar  to  the  one  in  Figure  A-l.B  can  be  used. 
Other  outputs  were  obtained  for  various  pole  locations. 

Results  obtained  for  the  following  cases  are  included,  with 
program  listing,  in  Appendix  A  and  B  respectively; 

1.  Complex  poles  with  moves  toward  a  constant 
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2.  Complex  poles  with  moves  toward  a  constant 

3.  Complex  poles  that  are  moving  toward  the  origin. 

4.  Real  negative  poles  with  one  of  which  moving  toward  the  origin 
while  the  other  remains  stationary. 

5.  Real  negative  poles  with  one  of  which  moving  toward  infinity 
while  the  other  remains  stationary. 

6.  Double  real  poles  moving  from  the  neighborhood  of  the  origin 
toward  infinity. 

7.  Pure  imaginary  poles  moving  toward  the  origin. 

8.  Pure  imaginary  poles  moving  away  from  the  origin. 

The  above  cases  are  now  described  in  the  following  sections.  In 
each  case  a  figure  is  included  showing  the  pole  locations  with  their 
values  shown  in  the  accompanying  table,  in  which  the  location  of  the 
corresponding  computer  output  plots  are  indicated. 
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Case  1 


This  case  illustrates  the  zeros  region  for  given  complex  poles  with 
a  moves  toward  a  constant  w  . 


Five  different  poles  (  conjugate  poles  are  not  shown  )  are  given  in 
accordance  to  show  the  zero  regions  for  each  as  tjp  remains  constant. 
The  locations  of  each  pole  are  tabulated  in  Table  1. 


Table  1 


Pole 

(13 

P 

Computer  Plot 

0.999 

1.0 

Figure  A-1.A,B 

^2 

0.707 

1.0 

Figure  A-2 

0.55 

1.0 

Figure  A-3 

0.25 

1.0 

Figure  A-4 

0.1 

1.0 

Figure  A-5 
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Case  2 


This  case  illustrates  the  zeros  region  for  given  complex  poles  with 
u)p  moves  toward  a  constant  • 

Five  different  pole  locations  are  shown  in  Table  2, 


Table  2 


Pole 

Computer  Plot 

C 

1 

1.0 

1.05 

All  plots  are 

^2 

1.0 

2.0 

shown  in 

C 

3 

1.0 

3.0 

figure  A-7 . 

C 

4 

1.0 

5.0 

1.0 

8.0 
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Case  3 


\ 


'  (  ju)) 

;\ 

S  ; 

S  ; : 

C4  I  ;  ,  ♦ 

'  1  1\ 

"^5  ;  :  i :  \ 

CJ3 

^Pa', 

''P4' 

^P4\ 

<^P5 

Figure  4. 


This  case  illustrates  the  zeros  region  for  given  complex  poles  that 
are  moving  toward  the  origin.  Five  different  pole  locations  are 
shown  in  Table  3. 


Table  3 


Pole 

T 

P 

0) 

P 

Computer  Plot 

4  .0 

30.652 

All  plots  are 

2.0 

15.326 

shown  in 

0.5 

3.834 

Figure  A-8. 

0.1 

0.766 

0.01 

0.077 

21 


Case  4 


This  case  illustrates  the  zeros  region  for  given  real  negative 
with  one  of  which  moving  toward  the  origin  while  the  other  remains 
t ionary. 


Figure  5. 


Five  different  pole  locations  are  shown  in  Table  4. 


Table  4 


Condition 

Pole 

Pole  P^ 

2 

8.0 

,  0 

4.0,  0 

6.0 

5.657 

s 

8  .0 

,  0 

1.0,  0 

4.5 

2.828 

S 

o 

00 

,  0 

0.5,  0 

4  .25 

2.0 

8  .0 

,  0 

0.1,  0 

4.05 

1.414 

8  .0 

,  0 

0.05,0 

4.025 

0.632 

All  plots  are  shown  in  computer  output  Figure  A-9. 


poles 

sta- 


Case  5 


This  case  illustrates  the  zeros  region  for  given  real  negative  poles 
with  one  of  which  moving  toward  infinity  while  the  other  remains  stati¬ 
onary. 


/fjtj 


P 


2 


<r 

p 


Figure  6 . 


Five  different  pole  locations  are  shown  in  Table  5 . 


Table  5 


Condition 

Pole 

Pole 

P 

2 

(T 

P 

‘^P 

C 

0 

0 

00 

■ 

2.0, 

0 

5.0 

4.0 

1 

C 

15.0,0 

2.0, 

0 

8.5 

5.477 

2 

C 

30.0,0 

2.0, 

0 

16.0 

•  7.746 

3 

"4 

50.0,0 

2.0, 

0 

26.0 

10.0 

C 

100.0,0 

2.0, 

0 

56.0 

14.142 

All  plots  are  shown  in  computer  output  Figure  A-10. 
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Figure  7. 


This  case  illustrates  the  zeros  region  for  given  double  real  poles 
moving  from  the  neighbourhood  of  the  origin  toward  infinity.  Five 
different  pole  locations  are  shown  in  Table  6. 


Table  6 


Condition 

Pole  P^-P2 

1 

u) 

P 

Computer  Plot 

1.0,  0 

1.0 

1.0 

All  plots  are 

^2 

2.0,  0 

2.0 

22.0 

shown  in  compu¬ 

5.0,  0 

5.0 

5.0 

ter  output 

10.0,  0 

10.0 

10.0 

25.0,  0 

25.0 

25.0 
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Case  7 


Figure  8. 

This  case  illustrates  the  zeros  region  for  given  pure  imaginary  poles 
moving  toward  the  origin.  Five  different  pole  locations  are  shown  in 
Table  7. 

Table  7 


Condition 

<r 

p 

cO 

P 

Computer  Plot 

C 

1 

0.0 

3.0 

All  plots  are 

<^2 

0.0 

1.0 

shown  in 

S 

0.0 

0.5 

Figure  A-12 . 

o 

o 

0.1 

0.0 

0.05 
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Case  8 


This  case  illustrates  the  zeros  region  for  given  pure  imaginary  poles 
moving  away  from  the  origin.  Five  different  pole  locations  are  shown  in 
Table  8. 


Table  8 


Condition 

a 

p 

Computer  Plot 

0.0 

4.0 

All  plots  are 

<=2 

0.0 

10.0 

shown  in  com¬ 

o 

o 

<.C .  0 

puter  output 

0.0 

100.0 

Figure  A- 13. 

0.0 

200.0 
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IV.  SENSITIVITY  STUDY 


A.  INTRODUCTION 

It  has  been  shown  [12],  that  in  a  linear  and  not  necessarily 
reciprocal  network,  with  only  one  independent  source  and  initially  at 
rest,  the  ratio  of  the  Laplace  transform  of  the  current  in  any  branch 
or  of  the  voltage  across  any  node  pair  to  the  Laplace  transform  of 
either  an  independent  voltage  or  current  source  is  given  by 
W  T  (0)  +  X  T  («) 

"W  =  --  w^V  - .  (1^) 


where 


X 

T(x) 


T  (0) 

X 


T  (») 

X 


an  adjustable  immitance. 

the  network  function  relating  a  response  to  an 
exitation. 

Lim  T(x)  =  the  network  function  evaluated  with 

X  — *■  o 

X  =  0,  is  assumed  to  be  finite. 

Lim  T(x)  =  the  network  function  evaluated  when 

X  — ►  ® 


X  =  ®,  is  assumed  to  be  finite. 

W  =  the  Thevenin  immitance  (independent  of  x)  seen  looking 
back  into  network  from  the  terminals  of  the  adjustable 
immitance.  W  is  assumed  to  be  finite. 

Equation  (16)  is  generic  in  form  in  that  x  can  be  either  an  impedance 
or  admittance.  In  this  paper  the  symbol  T(x)  is  used  to  denote  a  driving 
point  immitance. 
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B.  GRAPHICAL  REPRESENTATION  OF  T(x) 

Consider  equation  (16)  for  the  case  of  a  driving  point  immitance 
function  as  shown  in  Figure  10. 

Dividing  numerator  and  denominator  of  equation  (16)  by  W,  we  get 


T(x) 


X 

T  (0)  +  -  T  (oo) 

^  w  ^ 


X 

1  +  - 

W 


which  can  be  rewritten  as. 


T(x) 


X 

T^(oo)  (1  + - )  +  T^(0)  -  T^(co) 


X 

1  +  — 
W 


or 

T  (0)  -  T  (->) 

T(x)  =  T^(co)  ^ - 

X 

1  + - 

W 


(17) 


Let 


T(x)  =  T^(«)  +  V 


(18) 


where 


V 


T  (0)  -  T  («) 

X  X 

X 

1  +  — 

W 


(19) 
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T(x). 


Figure  10. 

Driving  p6int  immitance  with  variable 
component  x. 
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From  equation  (19)  we  have 

X 

V  +  V  -  »  T  (0)  -  T  (co)  (20) 

y  X  X 

we  note  that  ^nd  ^.re  fixed  quantities. 

Now  treat  each  quantity  in  (20)  as  phasor,  i.e. 

V  =  |v|  exp. 

and 

W  -  |w|  exp.  j62 

then  the  left  hand  side  of  equation  (20)  represents  the  sum  of  two  phasors. 

If  X  is  real,  then  the  phasor 
X  V 

V  -  =  X  I - 1  exp.  j(0  -  9  ) 

W  W  ^  ^ 

is  lagging  behind  the  phasor  V  by  an  angle 
9=8^-  9^. 

Since  fixed  quantities,  the  right  hand  side 

of  equation  (20)  is  constant  so  that  the  left  hand  side  of  (20)  is 
constant.  This  requires  that  the  angle  must  be  constant  (as  x  varies). 

This  implies  that  the  tip  of  phasor  V  describes  the  arc  of  the 
circle,  of  which  the  quantity  ^  chord  as  shown  in 

Figure  11. 

To  construct  the  circle  we  proceed  as  follows  [13] : 

1.  Plot  phasors  from  the  origin. 

2.  Draw  a  straight  line  through  the  tips  of  the  two  phasors  in 
step  1. 
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3.  Calculate  the  angle  and  draw  a  straight  line  from  the  tip  of 
T^(0)  forming  angle  with  having  the  value  -0  measuring  clock¬ 


wise  from  D. 


1 


Therefore  is  tangent  to  the  circular  locus  of  T^. 

The  centre  of  this  circle,  0*,  is  located  at  the  intersection  of 
the  two  perpendicular  lines  erected  on  the  chord  in  the  middle  of 
T^(0)  -  at  the  tip  of  shown  in  Figure  12. 

Knowing  the  locus  of  T(x)  as  x  varies  from  zero  to  infinity,  we 
can  graphically  determine  T(x)  at  any  value  of  x, 

C.  GRAPHICAL  REPRESENTATION  OF  SENSITIVITY  PHASOR 
If  we  define  the  sensitivity  as 

T(x)  ^  d  T(x) 

X  dx 

then  with  equation  (16)  we  can  obtain,  [13] 


gXCx) 


T^(0)  - 


^2 

X 


W  1  +  — 


W 


or 


V 


X 

V - 

W 


T^(0)  -  T^(«) 


where 


V 


Tx(0)  -  T^(«) 


X 


1  + - 

W 
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Figure  12- 


Graphical  representation  of  locus  of  T(x), 


Graphical  representation  of  sensitivity 
phasor  S, 


It  can  be  shown  [13]  that  the  sensitivity  of  a  function  may  be 
graphically  represented  by  a  phasor  such  as  PS  as  shown  in  Figure  13. 

D.  SAMPLE  ILLUSTRATION 

As  an  illustration  of  the  above  discussion  in  this  chapter,  let  us 
consider  a  biquadratic  p.r.  function  as  an  RC  driving  point  admittance, 
which  can  be  realized  in  Foster  2  form  as  shown  in  Figure  14. 

Let 


Y(s) 


2 

s  +  4s  +  3 

s^  +  6s  +  8 


(s  +  1)  (s  +  3) 

(s  +  2)  (s  +  4) 


which  can  be  realized  in  Foster  2  form  with  the  element  value  as  given 
in  Figure  5,  in  which  all  resistance  and  capacitances  are  measured  in 
Ohms  and  Farads  respectively. 

Now  let  R2  change  by  the  amount  AR,  such  that 

0  <  AR  <  w 

Thus,  Figure  15  is  modified  as  shown  in  Fgiure  16. 

We  can  calculate  T^(0)  ,  W  as  follows: 


T^(0)  =  Y(s,  AR=0) 


+ 


1 

+  1/sC^ 


+ 


1 

R^  +  l/sC^ 


07c  4.  «125  s  .  .09375  s 

.5  s  +  1  .25  s  +  1 


s^  +  4s  +  3 

s^  +  6s  +  8 


.  (24) 


T^(«)  =  Y(s,  AR=»)  =  —  +  r-Vi/Tc-- 

oil 


=  .375  + 


.125  s 

.5  s  +  1 


2.5  s  +  3 
4  s  +  8 


(25) 
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Figure  14. 

The  realization  of  Foster-2  form  of 
driving  point  admittance  Ks). 


Figure  15. 

The  realization  of  biquadratic  p.r.  driving  point 
admittance  in  Foster  2  form,  where 

Y(s)  =  _ 

s^+6s+8 
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w 


sC^  Y(s,  AR=co) 


.  .  10.667  .  4  s  +  8 

Z.OO/  +  - - -  +  -T—Z - —T 

I  Z.D  S  +  3 


14.5  s  +  40  s  +  99 

7.5  +  9  s 


(26) 


From  expression  (26)  we  obtain  the  expression  to  be 


X 

W 


"1.5  s  X  +  9  s  X 


(27) 


14.5  s  +  40  s  +  99 
Substitute  equations  (24),  (25),  (26)  and  (27)  into  equation  (22)  we  get. 


T(x) 


s  +  4s  +  3  2.5  s  +  3 

S^  +  6s  +  8_ 4  s  +  8 


14 


.5  +  40  s  +  99  ^  7.5  xs^  +  9  xs  ~|^ 

7.5  +  9  s  L  14.5  +  40  s  +  9?J 


,  (28) 


then  for  simplification  we  write  equation  (28)  above  as. 


T(x) 

S 

X 


N(s,x) 

D(s,x) 


(29) 


where  N9s,x)  and  D(s,x)  represent  the  numerator  and  denominator  of  the 
equation  (28)  respectively. 

In  the  frequency  domain,  the  expression  in  equation  (29)  can  be 
written  as , 


T(x) 

S  (jw,x)  = 

X 


N(.1cl),x) 

D(jw,x) 


(30) 
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Figure  16. 

Driving  point  admittance  Y(s)  with  variable  resis¬ 
tance  4R, 


1  w 

LdRi 

Li.. 

1 

1 

1 

1 

Y(s)  ;  ^ 

N.  R 

o 

1  i— 1 

J  i  1 

1 

1 

1 

1 

6 - 

1 

Figure  17. 

Driving  point  admittance  Y(s)  with  variable  resis¬ 
tance  R^'  =  ^2  ^ 
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I 

=  _  Re  lN(ja>,x)]  +  j  Im  [N(.ia),x)] 

Re  [D(j(D,x)]  +  j  Im  [D(ja>,x)] 

where  Re  [N(jw,x)],  Re  [D(ja),x)],  Im  [N(jw,x)]  and  Im  [D(j(ji),x)]  are 
real  and  imaginary  parts  of  equation  (30)  repectively. 

From  equation  (30a)  we  can  obtain  the  magnitude  of  that  sensitivity 
function.  If  we  write  equation  (30a)  in  exponential  form  we  have 


sT(x) 

X 


(jw,x) 


exp 


where 


sKx) 

X 


magnitude  of  equation  (30a) 


{Re  [N(ja),x)  ]  +  {lm[N(j(i),x)  ]  }^ 

(Re  [D(jco,x)]  ^  +  {Im[D(j(i),x)  ]  }^ 


(31) 


=  the  argument  of  the  numerator  of  equation  (30a)  . 

=  the  argument  of  the  demominator  of  equation  (30a) . 

E.  COMPUTER  STUDY  ON  SENSITIVITY 

In  this  section,  the  sensitivity  of  the  function  T(x)  with  respect  to 
the  variation  of  one  of  its  components  will  be  investigated  by  computer. 
Consider  the  general  expression  Y(x)  in  equation  (16),  which  is  repeated 
here  for  convenience. 

W  T  (s,x  =  0)  +  X  T  (s,  X  =  «>) 

T(x)  =  -  .  (16) 

W  +  s 


The  sensitivity  function  of  this  equation  is 


eT(x)  _ 


T^(s,  X  =  0) 


T  (s,  X  = 

X  * 


«  f  +  -f-] 


(22) 
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As  an  example  let  us  consider  the  function  Y(s),  as  we  derived  in  section  C 


Y(s) 


(s  +  DCs  +  3) 

(s  +  2)  (s  +  4) 


and  the  realization  of  this  function  in  Foster-2  form  is  shown  in  Figure 
17. 

For  our  investigation,  let  change  with  the  amount  of  AR  in  series 
with  R^  and  let 

R^'  *  ^2 

T  (s  R  * ) 

To  obtain  the  sensitivity  S  *  2  in  terms  of  frequency  and  R«’, 

^2  ^ 
we  use  the  same  calculation  from  equation  (30) . 

The  computer  output  in  Appendix  B  shows  us  the  plot  of  the  magnitude 

of  sensitivity  function  (30)  versus  ^2'*  this  purpose  we  use  equation 

(31)  as  reference. 

Various  cases  are  investigated  using  computer,  resulting  in  plots 
of  the  magnitude  of  equation  (31)  versus  R^’  for  different  frequencies. 
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V.  CONCLUSION 


In  this  paper,  a  comparative  study  has  been  made  on  the  realizability 
conditions  of  a  general  biquadratic  function.  The  graphical  method  of 
Chan  and  Phung  [10]  was  used  for  this  study  and  computer  programs  are 
written  for  verifications.  General  observations  can  be  made  as  follows: 

1.  The  realizability  for  the  zeros  expands  considerably  as  the  poles 
move  from  the  immaginary  axis  to  the  real  axis.  This  results  both  from 
C^  moving  away  from  the  origin  and  C^  moving  toward  the  origin. 

2.  The  limits  on  C^  and  C^  are  found  easily  from  equations  (3a)  and 

(3b) .  As  Op  approaches  zero,  corresponding  to  the  poles  moving  toward 

the  10)  axis,  both  C  and  C.  tend  toward  the  circle  o)  =  o)  for  all  (p  \ 
o  1  z  p  z 

As  can  be  seen  from  Figure  A-6,  C^  and  C^  separate  very  rapidly  for 

small  movements  of  the  poles  away  from  the  jm  axis.  We  can  conclude  that 

when  a  =  0  we  require  o)  =  o)  in  order  for  Z(s)  to  be  realizable  with 
P  ^  z  p 

complex  zeros.  This  can  also  be  seen  easily  from  equation  (9). 

3.  As  Op  approaches  u)^,  corresponding  to  the  poles  moving  toward 
the  real  axis,  C^  approaches  a  value  of  -5.828(Dp  on  the  real  axis,  and 

Cj^  similarly  approaches  a  value  of  -0.172a)p,  as  we  can  see  in  Figure  A-l.A. 

A  brief  study  was  made  on  sensitivity  functions  [13],  some  computer 
results  are  obtained.  However,  due  to  time  limitation,  a  complete 
investigation  in  sensitivity  was  not  possible  and  is  left  as  a 

suggested  topic  for  further  studies. 
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APPENDIX  A 


COMPUTER  OUTPUT 


A.  ZERO’S  REGION 

The  illustration  of  zeros  regions  will  be  shown  in  this  appendix. 

The  computer  outputs  show  us  the  regions  for  various  given  omega 
and  sigma  poles. 

Various  cases  are  explained  as  follows: 

Figure  A- la  and  A- lb  show  the  zeros  regions  for  given  sigma  pole 
=  0.999  and  omega  pole  =  1.0  using  Calcomp  and  PLOTP  in  the  computer 
program  respectively. 

All  the  following  plots  show  us  the  computer  outputs  by  using  Calcom 
in  the  computer  program. 

Figure  A-2,  A-3,  A-4  and  A-5  show  us  the  plots  of  zeros  regions  for 
given  sigma  poles  equal  to  0.707,  0.55,  0.25  and  0.1  respectively  with 
the  amount  of  omega  pole  equal  to  1.0  each. 

Figure  A-6  shows  all  plots  given  in  Figure  A-la,  A-2,  A-3,  A-4  and 
A-5  in  one  computer  output. 

Figure  A-7  shows  the  plots  of  zeros  regions  as  pole  moves  at  sigma 
pole  constant. 

Figure  A- 8  shows  the  zeros  regions  as  pole  moves  to  the  origin. 

Figure  A-9  and  A-10  show  the  plots  of  zeros  regions  as  pole  moves 
to  the  origin  and  to  infinity  respectively. 

Figure  A-11  shows  the  plot  of  zeros  regions  of  a  double  real  negative 
poles  as  both  move  along  the  negative  axis. 

Figure  A-12  and  A- 13  show  us  the  zeros  region  plots  of  a  given  pair 
of  pure  negative  poles,  as  both  move  to  the  origin  and  to  the  infinity 
respectively. 
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B.  SENSITIVITY  FUNCTION 


Various  cases  of  sensitivity  as  a  function  of  are  illustrated  in 
the  computer  outputs  for  different  particular  frequencies. 

The  computer  outputs  show  us  the  plots  of  sensitivity  versus  R*  ^ 
and  the  tabulations  of  this  plot. 

The  resistance  R’2  is  select  to  be  varied  from  the  amount  of  1  Ohm 
to  50  Ohms  with  increment  to  1  Ohm.  And  the  frequencies  are  choosen 
beginning  at  1  rad/sec.  up  to  10,000  rad/sec. 

Figure  B-1,  B-2,  .  B-9  show  us  the  sensitivity  functions  as 

R*2  moves  from  1  to  50  Ohms,  at  the  frequencies  of  1,  5,  10,  50,  100, 
500,  1000,  5000  and  10000  rad/sec.  The  tabulations  of  those  plots 
are  also  shown  in  Table  B-1,  B-2  . B-9  respectively. 

Figure  B-10  shows  us  five  different  plots  of  sensitivity  functions 
illustrated  in  one  computer  output. 
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The  plot  of  eensltivity  va  R^'  la  shown  for  given 
frequency  1  r-f’/aec. 
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Sensitivity  plot  vs  R^'  is  shown  here  for  given 
frequency  5  rad/sec. 
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for  given  frequency  10  rad/sec. 
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The  plot  of  eensltlvlty  ve  R^'  is  ehown  here  for  gi' 
frequency  50  rad/eec. 
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plot  of  eenBitlvity  vs  R^'  for  given  frequency 
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frequency  ^0  rad/eec. 
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frequency  1000  rad/sec, 
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The  plot  of  oeneltlvlty  vs  R. 
given  frequency  5000  rad/sec, 
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The  plot  of  Bensltlvit^  vs  R2' 
given  freouency  10000  rad/sec. 
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plots  of  sensitivity  vs  R. 
given  various  frequencies, 


The  V4LUE  CF  SENilTIVITT  i&T  GIVEN  W  = 


P2S 

SENSITIVITY 

C  .  0 

-0.005995 

l.GCO 

-0.005160 

2.000 

-0.004489 

2. COO 

-0.0C3940 

A. COO 

-0.003436 

5.0C0 

-0.002107 

6. COO 

-0.002736 

7.000 

-0.002512 

8.000 

-0.002277 

c.OOO 

-C. 002073 

1C. COO 

-0.001896 

11.000 

-0.001740 

12. COO 

-C.C01603 

12.000 

-0.001481 

lA.CCC 

-0.001373 

15.CCC 

-0.001276 

16. COO 

-0.001189 

17.000 

-0.001111 

18. COO 

-C. 001040 

1^.000 

-O.OOC976 

2C.CC0 

-0.000917 

21.000 

-C. 000864 

22.000 

-0.000815 

22.CC0 

-C.00C770 

2^. COO 

-C. 000729 

25.000 

-C.00C691 

26.CC0 

-C.000656 

27.000 

-C.000622 

28.000 

-0.000593 

2S.CCC 

-C.C00565 

2C.C00 

-0.000539 

21.000 

-0.000515 

22.CC0 

-C.GC0492 

22.000 

-0.000471 

2^.CC0 

-0.000451 

25. COO 

-C. 000432 

26.000 

-0.000415 

27.000 

-0.000398 

38.000 

-C.000382 

3S.CC0 

-0.000268 

AC. COO 

-C.000254 

Al.OOC 

-0.C00241 

A2.000 

-0 .000229 

A2.C0C 

-C. 000217 

AA.COO 

-0.000206 

A5.00C 

-0.000296 

A6.000 

-C. 000286 

47.000 

-0.000276 

48.000 

-0.000267 

4G.COO 

-0.000259 

5C.0C0 

-0.000251 

TABLE  B 
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THE  VALUE  CF  SENEITIVITV  AT  GIVEN  W 


5,00C 


R2S 

SENSITIVITY 

G.O 

-0.086060 

l.CCO 

-0.054744 

2.CC0 

-0.037862 

3,CCO 

-0.027736 

4.000 

-0.021188 

5.0CC 

-0.016  711 

6. COO 

-0.013516 

7.000 

-0.011156 

e.oco 

-0.009365 

G.COO 

-0.007972 

IC.QCO 

-0.006868 

ll.CCC 

-0.005979 

12.000 

-0.005252 

12.000 

-0.004649 

14.000 

-0.004145 

15.000 

-0.003718 

16.0CC 

-0.003354 

17. COO 

-0.003041 

18. COO 

-0.002770 

IS.CCC 

-C  .002534 

2C.0CC 

-C. 002326 

21.000 

-0.002143 

22.000 

-0.C01981 

23.000 

-C  .001837 

24.000 

-C. 001707 

25.000 

-0.001591 

26.000 

-C. 001487 

27.000 

-0.001392 

28.000 

-0.001306 

29.000 

-C. 001228 

30.000 

-0.001157 

31.000 

-C  .001091 

32.000 

-0.001031 

33.000 

-C. 000976 

34.000 

-C. 000925 

35.000 

-C.000878 

36.000 

-C.000835 

37.000 

-C. 000795 

38.000 

-0.000757 

39.000 

-0.000722 

40.000 

-C.00C690 

41.000 

-0.000659 

42.0CC 

-C. 000631 

43.000 

-C.C00604 

44.000 

-C. 000579 

45.000 

-0  .000556 

46.000 

-C.C00534 

47.000 

-C. 000513 

48.000 

-C.C0C494 

49.000 

-C.G00475 

5C.CCC 

-C. 000458 

TABLE  B  ->2. 
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THE  V^LUE  CF  SENSniVITT  /ST  GIVEN  W 


10.000 


P2S 

SENSITI VITY 

0  •  c 

-C. 227638 

l.CQO 

-'C.122C24 

2.0C0 

-0.076772 

2.CC0 

-0.052258 

^•CQO 

-0.037941 

5.CQQ 

-0.028735 

6. COO 

-0.022506 

7. COO 

-0.018097 

8.0CQ 

-0.014865 

«.CQC 

-0.012425 

IC.OCO 

-0.010540 

ll.COQ 

-0.009052 

12.0CC 

-C.0C7858 

13.000 

-0.006885 

lA.CCC 

-C. 006082 

15.000 

-0.005412 

16. COO 

-C  .004846 

17.000 

-C. 004365 

18.0C0 

-C. 002952 

IS. 000 

-0.003595 

20.CCC 

-0.003284 

21.000 

-0.003012 

22.000 

-0.002772 

22.000 

-0.002560 

2^.000 

-C. 002271 

25.000 

-0.C02202 

26.CC0 

-C.002C51 

27.000 

-0.001915 

28.000 

-0.001792 

2S.000 

-0.0C1680 

20. COO 

-0.001579 

21.000 

-0.001486 

22.CC0 

-O.C01401 

23.000 

-C  .001224 

24.000 

-0.001252 

25.0C0 

-0.001187 

26. COO 

-0.001126 

27.0C0 

-0.001070 

28.CC0 

-0.001018 

2S.CC0 

-C. 000969 

4C.0Q0 

-0.000924 

41.000 

-0.000882 

42.0C0 

-0.000843 

42. COO 

-C.C00807 

44.000 

-0.000772 

45.CC0 

-C.C0C740 

46.CC0 

-C.CCG71Q 

47.000 

-C  .000682 

48. COO 

-0.000655 

4S.CC0 

-0.000620 

5C.CCC 

-0.000606 

TABLE  B  -  3 
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THH  VALUE  CF  SENSITIVITY  AT  GIVEN  W  =  50.000 


P2S 

SENSITIVITY 

c  •  c 

-C.67A75a 

1.000 

-0.270A19 

2.000 

-0.1AA071 

3.000 

-0.C89189 

4.CC0 

-C. 060562 

5.000 

-C.0A3783 

6. COO 

-0.032117 

7.000 

-C.025S19 

8.CC0 

-0.020626 

S.COO 

-0.017112 

IC.OCO 

-0.01A3C5 

ll.CCO 

-0.012135 

12.000 

-C.010A2A 

13.000 

-0.009050 

14.CC0 

-0.007932 

15.0CC 

-O.CCTCCa 

16. COO 

-0.006237 

17.000 

-0.005587 

le.coo 

-C.C05023 

IS. 000 

-0.00A557 

2C.CC0 

-0.00A1A6 

21.000 

-0.003788 

22.CC0 

-0.003A75 

23.000 

-0.003199 

24.CC0 

-0  .C0295A 

25. COG 

-C. 002737 

26. COO 

-0.0025A2 

27.000 

-0.002368 

ze.ccc 

-C. 002211 

2S.0C0 

-C. 002069 

2C.0C0 

-0.0019A1 

31.CC0 

-C.C0ie2A 

32. COO 

-C. 001717 

33.000 

-0.001619 

2A.CC0 

-0.001530 

35.000 

-0.001AA7 

36.CC0 

-0.001371 

37.CC0 

-0.001301 

36.CCC 

-C. 001227 

2S.CCC 

-0.001176 

AO. COO 

-0.001121 

Al.OCC 

-0.001069 

A  2 . 0  CO 

-0.001020 

A3.CC0 

-C. 000975 

AA.OOO 

-0.000922 

A5.C00 

-C.C00893 

A6.000 

-0.000856 

A7.0C0 

-0.000821 

A8.000 

-0.000789 

AS.CCC 

-0.000758 

5C.CCC 

-0.000729 

TABLE  B  -  4 
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THE  VALUE  CF  SENSITIVIIV  AT  GIVE^  W 


100,000 


RZS 

SENSITIVITY 

0,0 

-0.739854 

1.000 

-0.284703 

2,0CC 

-0.149352 

3,000 

-0.091733 

A. COO 

-0.061999 

5.00C 

-0.044683 

6. COO 

-0,033724 

7.CC0 

-0.026352 

8.CC0 

-0.021157 

c.CCO 

-0,017359 

10.000 

-0.014499 

11.000 

-0.012292 

12.CC0 

-C.01C553 

13,000 

-0.009158 

l^.CCO 

-C.C08C23 

15.000 

-0.C07C66 

16.000 

-0,006304 

17.000 

-0,005645 

lE.QCO 

-0.005084 

l^.OCC 

-0.G04603 

2C.C00 

-0.004187 

21.000 

-0.003625 

22.000 

-0.003508 

23.000 

-C. 003229 

24 .000 

-0,002981 

25.000 

-0,002762 

26.000 

-0.002565 

27.000 

-0.002389 

26.CCC 

-0,002230 

2S.C00 

-C.C02C87 

30.000 

-0.001957 

31.000 

-C.G01839 

32.000 

-0.001731 

33.000 

-0.001633 

34.000 

-0,001542 

35.000 

-0.001459 

36.000 

-C.C01383 

37.000 

-C. 001312 

38.000 

-0.001247 

3^.000 

-C. 001186 

40.000 

-C.COl 130 

41.000 

-0.001077 

42,000 

-0.001028 

43. COO 

-C.0C0983 

44,000 

-O.OC0940 

45.000 

-C.CC0900 

46.000 

-0  .000863 

47.000 

-C. 000828 

48.CCC 

-C.C0C795 

4C.0C0 

-0.00C764 

50. COO 

-C. 000734 

TABLE  B  ~  3 
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THE  VALLE  CF  SENSITIVITV  AT  GIVEN  M  =  500.000 


R2S 

SENSITIVITY 

0.0 

-0.764546 

l.CCO 

-0.289726 

2.0CO 

-0.151158 

2. GOO 

-0.092591 

4. COO 

-0. 062480 

5.0C0 

-0.044982 

6. COO 

-0.032925 

7. COO 

-0.026495 

8.000 

-0.021263 

^.CCC 

-0.017441 

IC.OCC 

-0.014563 

ll.CCO 

-0.012244 

12.000 

-0.010595 

13.CCC 

-0.C09193 

14.0C0 

-0.008053 

15.CCC 

-0.007112 

16.C0C 

-0.006227 

17.0C0 

-0.005665 

16.000 

-0.005101 

19.CC0 

-0.004618 

2C.CC0 

-0.004200 

21.0CC 

-0.003637 

22. COO 

-0.003519 

23.000 

-0.002238 

24.000 

-C.0C2990 

25.CC0 

-O.C02770 

26.000 

-0.002573 

27.000 

-0.002296 

26.000 

-0. 002237 

2^.000 

-0.002093 

50. COO 

-0.001963 

31.000 

-0.001644 

32.000 

-0.001736 

22.000 

-0.001637 

24.000 

-0.001546 

35.000 

-0.001463 

26.000 

-0.001286 

27.000 

-0.001215 

38.000 

-0.001250 

3S.C00 

-0.001189 

40. COO 

-0.001132 

41.000 

-0.001080 

42.000 

-0.001031 

42.000 

-0.000985 

44.000 

-0.000942 

45.000 

-0.000902 

46.000 

-C.OOG865 

47.000 

-C  .000620 

46.000 

-0.000797 

49.000 

-0.000765 

50.000 

-0.000726 

TABLE  B  -  6 
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THE  VALUE  CF  SENSITIVITY  AT  GIVEN  W  =  1000.000 


P2S 

SENSITIVITY 

€.0 

-0.765354 

l.OCO 

-0.289887 

2.CC0 

-0.151215 

3.CC0 

-0.092619 

^•COO 

-0.062495 

5.C0C 

-0.044992 

6.000 

-0.033921 

7.QC0 

-0.026499 

6.CCC 

-0.021266 

S.OCO 

-0.017^43 

1C. 000 

-0.014565 

11. COO 

-C.C12345 

12.000 

-0.010596 

13.CQC 

-0.009195 

l^.OCO 

-C.008C54 

15.0CC 

-C.C07113 

16.0C0 

-0.006327 

17. COO 

-0.005665 

18.000 

-0.005102 

IS.CCO 

-0.004619 

2C.CCC 

-C.CG4201 

21.000 

-0.003827 

22.CCC 

-0.003519 

23.CC0 

-0.003239 

2^.CCC 

-0.002991 

25.CCC 

-C.CC2770 

26.CCC 

-0.002573 

27.000 

-0.002396 

28.000 

-0.002237 

2G.C00 

-0.002093 

3C.CCC 

-0.001963 

31.000 

-0.001844 

32.000 

-0.001736 

33.CCC 

-0.001637 

34.C0C 

-0.001546 

55.CC0 

-0.001463 

36. COO 

-0.001386 

37.000 

-0.001315 

38.000 

-0.001250 

39.000 

-0.001189 

40.000 

-0.001133 

41.000 

-0.001080 

42.000 

-0.001031 

43.000 

-C. 000985 

44.000 

-0.000943 

45. COO 

-0.000903 

46.000 

-0.CO0865 

47.000 

-C. 000830 

48.000 

-0.0CC797 

49.000 

-0.000765 

5C.CCC 

-0.000736 

TABLE  B  >  7 
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THE  VALUE  CF  SENSITIVUV  AT  GIV£^  W  =  5000.000 


P2S 

SENSITIVITY 

c.c 

-0.765613 

l.COC 

-0.289929 

2.0GC 

-0.151224 

2.00C 

-0.092627 

4.CC0 

-C. 062500 

5.CCC 

-0.044995 

6.CCC 

-0.033923 

7.CC0 

-C. 026501 

8.000 

-0.021267 

9.QC0 

-0.017444 

IC.CCC 

-0.014566 

11.000 

-0.012246 

12.0CC 

-0.010597 

12.000 

-C. 009  195 

14.0CC 

-0.008054 

15.000 

-0.007113 

16.CCC 

-0.C06327 

17.0CC 

-0.005665 

16.000 

-0.005102 

Ic.CCO 

-0.004619 

2C.CCC 

-0.004201 

21.CC0 

-0. 003627 

22.000 

-0.002519 

23.000 

-0.002239 

24.000 

-0.002991 

25.000 

-0.002770 

26.CCC 

-0.002573 

27.000 

-0.002396 

28.000 

-0.002237 

29.000 

-C. 002093 

20. COO 

-0.001963 

21.CC0 

-0.001844 

22.000 

-0.001736 

33.000 

-0.001637 

24.000 

-0.001547 

25.000 

-0.001463 

26.000 

-0.001286 

37.000 

-0.001315 

26.CC0 

-0.001250 

39.000 

-0.001189 

40.000 

-0.001133 

41.000 

-c.ocioao 

42.000 

-0.001021 

42.000 

-0.000985 

44.000 

-0. 000942 

45.000 

-C. 000903 

46.000 

-0.000665 

47.000 

-C.CC083C 

46.000 

-0.C0C797 

49.000 

-0.000766 

50.000 

-0.000736 

TABLE  B  -  8 
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THE  VALUE  CF  SENSITIVm  AT  GIVEN  W  =  10000.000 


R2S 

SENSITIVITY 

C  •  0 

-0.765621 

I.OOO 

-0.289940 

2.0CC 

-0.151234 

3.CC0 

-0.092627 

4.0C0 

-0.062500 

5.CC0 

-0.044995 

6.CG0 

-0.033933 

T.OCO 

-0.026501 

8.000 

-0.C21267 

<:.CCO 

-0.017444 

IC.CCC 

-0.014566 

11.000 

-0.012346 

12.000 

-C. 010597 

13.CC0 

-0.009195 

14.000 

-0.008054 

15. COG 

-O.OC7113 

16.CCC 

-0.006327 

17.0C0 

-0.C05665 

18.CC0 

-O.C05102 

l^.CCO 

-0.CC4619 

2C.CCC 

-C.C04201 

21.000 

-C .003837 

22.000 

-0.003519 

23. COO 

-0.C03239 

24.CCC 

-0.002991 

25.000 

-C. 002770 

26. COO 

-0.002573 

27. COG 

-0.002396 

28. COG 

-C.C02237 

29. COO 

-0.002093 

3C.CCC 

-0.001963 

31.000 

-  -C.CC1844 

32.0CC 

-0.001736 

33.CC0 

-0.001637 

34.000 

-0.001547 

35.000 

-0.001463 

36.CC0 

-0.001386 

37.CC0 

-C.C01315 

38.0C0 

-0.001250 

39.000 

-0.001189 

4C.CC0 

-0.001133 

41. COO 

-C. 001080 

42.CCC 

-0.001031 

43.000 

-0.000985 

44.CQC 

-C.0C0943 

45. COO 

-O.C00903 

46. COO 

-0  .OCO  86  5 

47. COO 

-C. 000830 

48.000 

-C. 000797 

49. COO 

-0.000766 

50. COO 

-C.OOC736 

TABLE  B  -  9 
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APPENDIX  B 


COMPUTER  PROGRAM 


A.  ZEROS  REGION 

This  program  will  plot  the  realizability  region  of  general  bi¬ 
quadratic  function  of  the  form 

+  2*SIGMAZ*s  +  OMEGAZ^ 

Z(s)  *  - 

+  2*SIGMAP*s  +  OMEGAP^ 

The  input  to  this  program  is  via  punched  card.  One  card  is  required 
for  each  plot.  The  input  variables  are  Sigmap  and  Omegap  specified  in 
columns  1  -  20  as  two  fields  of  F  10.3  format. 

The  output  is  a  plot  of  allowable  region  for  zeros  of  Z(s)  for  a  given 
pair  of  poles.  The  criterion  for  the  plot  is  that  Z(s)  must  be  positive 
real. 

The  output  may  be  either  one  (or  both)  of  two  types  of  plots  depend¬ 
ing  on  the  needs  of  the  programmers.  If  a  fast  solution  is  desired  the 
"CALL  DRAW"  statements  must  be  removed  to  enable  the  program  to  be  run 
on  QUICKRUN.  The  result  will  be  a  rough  plot  on  the  line  printer 
output. 

If  TnayTmiim  accuracy  is  desired  a  CALCOMP  plot  may  be  obtained  by 
making  a  regular  FORTCLG  run  with  the  "CALL  DRAW"  cards  in  the  deck. 

In  this  case  a  line  printer  plot  will  also  be  generated,  unless  the 
"CALL  PLOTP"  cards  are  moved-  For  the  CALCOMP  run,  "REGION. GO— 10 OK" 
must  be  specified  on  the  EXEC  card. 

The  routines  DRAW  and  PLOTP  are  called  by  this  program.  Both  are 
in  the  MPSLIB  at  NPS.  Neither  a  copy  of  the  decks  nor  special  job 
control  language  is  needed  to  run  this  program  on  the  NPS  IBM  360. 
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Computer  program  for  the  realizability  region 
of  general  biquadratic  function. 


The  computer  program  is  shown  in  Figure  A- 14. 


B.  SENSITIVITY  FUNCTION 

This  program  will  plot  the  sensitivity  of  general  biquadratic 
function  with  respect  to  the  changing  of  its  components  at  a  particular 
frequency. 

The  equation  used  in  this  program  is. 


(joi.Rp 


W  (T(ja),R^  -  0)  -  T(ju,R^  =  INF.) 

(W  +  R^)^ 


where  T(ja),Rp,  W,  Kju.R^  =  0),  T(ja),R2  =  INF.)  and  R^  correspond  to 
Figure  17  in  Chapter  IV.  E. 

The  output  is  a  plot  of  the  magnitude  of  the  sensitivity  function 
above  versus  R^* 

The  input  of  this  program  is  via  punched  card.  One  card  is  required 
for  each  plot.  The  input  variable  is  frequency  in  rad/sec.  specified 
in  columns  1  -  10  as  one  field  of  F  10 .*3  format. 

The  output  is  a  rough  plot  on  line  printer  output. 

The  routine  "CALL  PLOT?"  is  called  by  this  program. 

The  computer  program  is  shown  in  Figure  B-11. 
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Computer  program  for  sensitivity  function  of 
given  particular  biquadratic  function. 
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